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A negative energy electron will have less energy the 
faster it moves and will have to absorb energy in 
order to be brought to rest. No particles of this nature 
have ever been observed. - - P. A. M. Dirac, 1930 



Abstract 

We study the longitudinal instabilities of two interpenetrating fluids interacting only 
through gravity. When one of the constituents is of relatively low density, it is possible 
to have a band of unstable wave numbers well separated from those involved in the 
usual Jeans instability. If the initial streaming is large enough, and there is no linear 
instability, the indefinite sign of the free energy has the possible consequence of explosive 
interactions between positive and negative energy modes in the nonlinear regime. The 
effect of dissipation on the negative energy modes is also examined. 
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1 BACKGROUND 



It used to be thought that the stellar component of two spiral galaxies would pass right 
through each other in the event of a collision and that only the gaseous components would 
merge. However, simulations over the past twenty years or so 1 ' 2 have shown that the macro- 
scopic energy of such a collision is quickly converted into internal energy and that merger 
of the stellar systems is a common natural outcome of a collision. How is this conversion 
effected? 

An answer to this may lie in the physics of streaming instabilities. In the context of 
plasma physics, interpenetrating electron beams produce the two-stream instability 3 whose 
gravitational analog has long been recognized, beginning with the investigations of Sweet 4 and 
Lynden-Bell. 5 Of course, in the case of galaxy collisions, which occur quickly, the conventional 
two-stream instability may operate too slowly to be effective. However, it is known that even 
when two streams of interacting particles are not linearly unstable, they may collectively 
produce negative energy modes that lead to an explosive nonlinear growth of perturbations 
for arbitrarily small disturbances. There are well-developed criteria for the occurrence of 
this explosive growth in plasma physics 6 and, as Lovelace et al. 7 have suggested in their 
analysis of counter-rotating galaxies, we may expect something analogous in the gravitational 
setting. Our aim here is to briefly develop this topic of negative energy modes for the case 
of gravitational interaction in the expectation that the phenomena involved will be found 
significant in a variety of astrophysical processes. 

Even in the event of linear instability, the case of counterstreaming populations is signif- 
icantly different from the standard gravitational instability, which occurs only for perturba- 
tion scales greater than the Jeans length. 8 When there are two interpenetrating fluids such 
as stars and gas, modes of arbitrary wavelength can be rendered unstable. Numerous inves- 
tigators have reported on these issues. Most of them (such as Ikeuchi et al., 9 Fridman & 
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Polyachenko, 10 and Araki 11 ) focused primarily on the symmetric situation of identical fluids 
in counterstreaming motion. In that case one finds that the spectrum of any instabilities aris- 
ing from the relative motion is wholly contained within the Jeans instability band, and this 
blurs the distinction between the two processes. This need not be true when this symmetry 
is broken, and indeed not all authors restricted themselves entirely to the symmetric case. 
The present venture into the asymmetric problem is intended to focus on the possibility of 
well-separated instability bands, which has not been elucidated in the gravitational context, 
as far as we are aware. 



To see the problem in its simplest version, it is useful to have a uniform medium as the 
unperturbed state. Rather than formulate the problem inconsistently to achieve this end, as 
Jeans did, we prefer the Einstein device of introducing a cosmological repulsion term. In the 
Newtonian setting we readily see how to redefine the gravitational potential so that, instead 
of introducing such a repulsion term, we fill space with a fluid of negative gravitational mass 
of density p A . As in the one-fluid plasma model, we treat this density as a constant since its 
purpose is to allow a gravitationally neutral background state. One may also contemplate 
the analog of the two-fluid plasma model in which this background antigravitational fluid has 
its own dynamics, but we do not do that here. The two dynamically active fluids we consider 
are gravitationally ordinary and polytropic. Thus, the Poisson equation is written here as, 



where V is the gravitational potential, p\ and pi are the source densities of the conventional 
fluids, and — p A is the cosmological background density. 
The equations of motion for the two fluids (j = 1,2) are, 



2 EQUATIONS OF MOTION 



W = 4ttG ( Pl + P2 - p A ) , 



(1) 



-Vpj - pjVV 



(2) 
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dtPj + V-ipjUj) = 0, (3) 

where we do not sum over repeated indices. Each fluid has a sound speed, c| =^~, and 
a Jeans wavenumber, kj- —^E^Ei^ where the hat signifies the equilibrium value and an 
uncustomary factor of 2 appears in the definition of the Jeans wavenumbers. 

We shall use natural units with /c" 1 as the length scale and (k J1 ci)^ 1 as the time scale. We 
further simplify the description by considering only longitudinal motions in one-dimension, so 
each single-component velocity field may be expressed as the gradient of a velocity potential: 
Uj = d x (j)j. The fundamental equations (l)-(3) take the dimensionless form, 

M 1 d t <h. + lM^d x <l> 1 ) 2 + ^— + V = B 1 (4) 

cM 2 d t 4> 2 + \c 2 Ml(d x 4>2? + + V = B 2 (5) 

dtpi + M&ipxdM = (6) 

dtP2 + cM 2 d x {p 2 d x (j>2) = (7) 

^-^-^P 2 + f = 0, (8) 



where the Mj =^ are Mach numbers, Uj measures the initial streaming velocities, f3 
c =— , and the Bernoulli constants B. are chosen to balance the basic state. 
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3 LINEAR THEORY 

We perturb from the state of uniform densities and constant velocities by setting (f>j = 
(— iy +1 x + 5(f) j, pj = 1 + Spj, and V = V + SV. The density terms of the Poisson equation 
(8) combine to vanish and V is a constant. Since the linearized equations are separable we 
may decompose the perturbations into normal modes proportional to exp(iwi — ikx) to find 
the dispersion relation, 

r( "' k) = 1 + 2 [(„ - ku.y rpj + 2 [ ( c + c^m 2 )^ - cW] ' 1 + ri + r2 = °' (9) 
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The quantity T, which we call the diagravic function by analogy with the dielectric function 
of electrodynamics, measures the collective response of the fluid to disturbances in the grav- 
itational field and will serve to indicate the energy signature of any normal mode (section 
4). 

For real k the solutions of (9) with complex oj correspond to instability; if w is real, then 
solutions of (9) with complex k can give rise to wave amplification instability. Here we analyze 
only the case of real k. However, we have to deal with both the traditional Jeans instability 
as well as the two-stream instability, the latter of which involves a sympathetic bunching 
of particles and is effective for creating instability when the phase speed of the disturbance 
conspires to create a resonance between different modes. 

3.1 Symmetric Case 

If both fluids have the same basic properties (c = 1 and [3 = 1), a frame exists in which 
M 1 = M 2 = M. The dispersion relation then simplifies into a manageable biquadratic with 
solutions, 

<4 = ~ + k 2 {M 2 + 1) ± yjj - 2k 2 M 2 + 4M 2 *; 4 . (10) 

For M = 0, we recover a simple version of the previously studied two-fluid Jeans 
problem. 10,12 ' 13 We find u\ = k 2 , corresponding to sound waves at all k, and lj 2 , = k 2 — 1, 
which is the conventional Jeans dispersion relation. The new acoustic modes arise because 
the aggregate fluid now allows motions unaffected by the gravitational field; for these modes 
the perturbed gravitational potential is zero. 

With relative velocity in the subsonic regime (0 < M < 1), there is only a single unstable 
mode that branches continuously from the Jeans mode at M = 0. This mode is unstable 
for all wavenumbers below a critical value that approaches infinity as M tends to unity from 
below (see Table 1). To study this limit, we let M 2 = 1 — a/k 2 with < a < 1. As 
k — > oo we find the approximate solution u; 2 ~ — — ^, which reveals a weak instability 
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at large k. Thus, weak relative streaming allows gravitational instability at arbitrarily small 
wavelengths. These large-/c instabilities do not arise for Maxwellian velocity distributions 
within the context of the Vlasov equation. 9 

For supersonic motion (M > 1) the large- A; gravitational instability is no longer present, 
but a new instability emerges that we call a two-stream instability since it owes its presence 
to the energy contained in the initial streaming motion. As M ranges from 1 to oo, the 
critical wavenumber for instability increases from k cr n = 1/2 to k cr n = \/2/2. 

The upper half of figure 1 shows that near k = the two-stream modes are wholly 
contained within the Jeans band. This fact coupled with the larger growth rates of the 
Jeans modes has led some to believe that the two-stream instability is swamped by the Jeans 
instability and is essentially unimportant. 11 As k increases the two-stream and Jeans modes 
collapse upon each other and together bifurcate into growing and damped oscillations. At 
still larger wavenumbers all motions are stable, propagating waves. The critical wavenumbers 
below which growth is possible at any Mach number are shown in Table 1. 

TABLE l 1 



Mach Range 


Mode Type 


k 2 

^crit 


limM^i k' 2 crit 


ImiAf^oo k 2 rit 


< M < 1 


Jeans 


1 

1-M 2 


oo 


Not Applicable 


1 < M 


Two-Stream 


VM 2 -1 


1 





4¥(M 2 -1+MVM 2 -1) 


4 


1 < M 


Jeans 


VM 2 -1 


1 


l 


4M(1-Af 2 +MVM 2 -1) 


4 


2 



3.2 Asymmetric Case 

When we relax the constraint of identical conditions in the two fluids, one of the more 

interesting consequences is the possibility of large wavenumber two-stream instability bands 

lr The value of k\ rit for the two-stream modes with M > 1 more accurately refers to the k value at which 
the unstable two-stream and Jeans branches merge. 



6 



well-separated from the Jeans instabilities clustered at small k. For illustration we consider 
the effect of changing the initial relative streaming Mi + cM.2 for fixed (3 and c. In fact, (3 
turns out to be the crucial parameter in achieving the spectral separation; variations in the 
sound speed ratio c widens both bands together. The distancing of a bubble of two-stream 
modes from the Jeans band is illustrated in figure 2. 

The large-fc two-stream modes can be explained qualitatively by examining the separate 
pieces of the dispersion relation (9). Since the streaming instability is related to resonant 
motions, it is revealing to examine the solutions to 1 + Y\{uj, k) = 1 + T2(u}, k) = in isolation 
and see where the curves intersect. The frequencies of these non-interacting modes are given 
by, 



-•j - A-Mi^-I (11) 



lo 2 = -ckM 2 ± ] jc 2 k 2 -| (12) 

When oj\ = UJ2, the assumed independent frequencies match one another for wavenumbers 
satisfying, 

k 2 -\ + Jc 2 k 2 - § 




2 V 2 

In the symmetric case where Mi = M2 = M, c = (3 = 1, we see that two of the resonances 
are lost except in the irrelevant cases M = and k = 0. The other pair of possibilities, 



kM = ±\Jk 2 — \, just restate the critical wavenumber condition for what we know to be the 
modified Jeans instability when M < 1. In the general case, we can expect another pair of 
intersections that account for the two-stream bubbles of figure 2. 
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4 NONLINEAR THEORY 



4.1 Hamiltonian Formulation and Energy Signature 

The dynamical equations (4)- (8) derive from a variational principle and a conserved 
Hamiltonian functional. The variational formulation has the advantage of shedding light on 
the relation between the energy content of the disturbances and nonlinear stability. Here 
we present results for the symmetric case, though the formalism follows through for the 
asymmetric case as well. 

The Hamiltonian and associated equations are, 

1 2 f L ( M 2 \ 
H = mT,J q dx^—pj^ + pjUj-VZj, (14) 

dth = {<f>j,H} = -— (15) 
()H 

dt Pj = { Pj ,H} = —, (16) 

7-1 

P 

where Uj(pj) = ^^_^ is the internal energy for the j th fluid and the Poisson bracket is 



defined by, 



For definiteness we have chosen a box-geometry of length L. 

The energy content of a particular mode when the perturbation amplitude is small is 
given by the second variation of H evaluated at equilibrium (the free energy). After some 
calculation this is seen to be, 

S 2 H = ^J q L dx (M 2 [5(f> 2 lx + 54> 2 X ] + 2M 2 [5 Pl 5<p lx - Sp 2 S ( (> 2x } + 5 P 2 + bp 2 - 25V 2 ) . 

(18) 

It may be verified that this functional is conserved by the equations of motion. Suppose we 
now insert into (18) an eigenfunction corresponding to a single stable mode with Im(a>)= 0. 
Employing overbars to denote eigenvector components and * for complex conjugation, we 



write, 

SV = Ve i{u}t - kx ^ + F* e -^**- fca: \ (19) 

and similarly for the other perturbation variables. Upon effecting the integrations, we can 
make use of the dispersion relation, 

r(w ' k) = 1+ 2[(u-kM)*-k*] + 2[{u + kMY-^\ = ° ' (20) 
to express the modal free energy in the compact form, 

5 2 H = -2Lk 2 | V | 2 u^. (21) 
duo 

Wherever 0Jj^-< a positive energy mode (PEM) is implied by (21), while the condition 



cj|^> defines a negative energy mode (NEM). This possibility of modes of either signa- 



ture has been elucidated in the plasma physics literature 14 and a gravitational analog was 
suggested by Lovelace et al. 7 in the context of thin, counter-rotating stellar disks. 

Figure 3 shows the diagravic function for both a subsonic and a supersonic case of stable 
modes with k = 3. In the subsonic regime, we find the Hamiltonian to be positive definite 
near equilibrium since ^f^< at every crossing of T on the w-axis. Right at the border 
of supersonic streaming (M = 1), concomitantly with the appearance of the two-stream 
instability, the T curves undergo a topological transition that allows the coexistence of positive 
and negative energy modes. The NEMs are slow modes in that they have smaller frequencies 
than their PEM counterparts. This is the typical situation; u> must pass through zero if 
the energy signature changes. 15 We expect from the precedents of plasma physics 6 that the 
simultaneous presence of positive and negative energy modes has dramatic consequences on 
the nonlinear stability of the system. 

4.2 Reduction to Action Angle Variables 

In the rest of this section we will concentrate on the nonlinear interactions between 
linearly stable modes in the symmetric problem (see figure 3). From a physical standpoint, 



attention is focused on situations where the disturbances are of sufficiently small scale so that 
the Jeans instability can be ignored, though there are no compelling reasons why this ought 
to be the case. We will further assume supersonic motion in order to examine the interaction 
of positive and negative energy modes, a situation we expect to be the most interesting. 
Under these assumptions the equations of motion achieve their simplest form in action-angle 
coordinates that we now develop. 

First we Fourier transform the field variables: 

Pi= E P m \t)e ik -* , - k m = 2 -^. (22) 

ra=— oo 

We can then write the free energy in terms of real variables as, 

* 2H =2 E (q T Aq + P T Bp), (23) 
m=l 

where q = (q\, q 2 , <73, q&) T , p = (pi,P2,P3,P4) T ', are linear combinations of the complex modal 
amplitudes, A and B are symmetric matrices (given in Casti 16 ), and the "T" indicates 
transpose. 

Defining the configuration variables z = (qi, . . . ,Pi) T , we recast the linearized equations 
in the form, 

|^v^ sLz . l^(_° a J) , a = (V„). <"> 

where J is the canonical 8x8 cosymplectic form. The next order of business is to construct 
a symplectic transformation that puts 5 2 H in its normal form. 17 ' 18 This can be achieved by 
writing z = SZ, where the matrix S consists of suitably ordered eigenvectors of L satisfying 
the symplectic condition, S T JS = J. After a final transformation to action-angle coordinates, 

the free energy expression (23) becomes a superposition of harmonic oscillators, 

00 

S 2 H = E (u+J\ + uj + J 2 uj^J a ) . (25) 

m=l 

The free energy is thus manifestly composed of two pairs each of positive and negative energy 
modes. 
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4.3 Three- Wave Resonance and Explosive Growth 

Energy conservation forbids nonlinear runaway growth if H is definite (Dirichlet's the- 
orem), as is the case here for subsonic motion. When the relative streaming is supersonic, 
however, interacting PEMs and NEMs can circumvent this restriction since they contribute 
energy of opposite sign. 

We demonstrate the possibility of explosive growth with a three-wave resonant interaction 
between two NEMs and one PEM. Since the energy signature of a mode is not Galilean 
invariant, the existence of a reference frame in which all three modes have the same signature 
implies nonlinear stability. It may be shown that there is no reference frame in which all three 
modes have the same energy signature if the highest frequency wave has opposite signature 
to that of the other two. 19 This provides a criterion for three- wave interactions leading to 
instability. 

The third-order resonance conditions for a triplet of modes are, 
m 1 k 1 + m 2 k 2 + m 3 k 3 = 

m 1 uj 1 ± m 2 uj 2 ± m 3 u 3 = (26) 
| m 1 | + | m 2 | + | m 3 | = 3 (m 1 ,m 2 ,m 3 integers). 

which here may be satisfied by (m 1 ,m 2 ,m 3 ) = (1,1,-1), (k 1 , k 2 , k 3 ) = (k m ,k m ,2k m ), and 
(uj 1 ,uj 2 ,lu 3 ) = (lj+,u)-,u-), where the u j are taken to be positive. One can see from figure 
1 that lo 1 > u 2 ,co 3 , so the relative signatures of this triplet are immune to a Galilean shift. 
Note from figure 1 that a resonant triplet involving two PEMs and one NEM would not have 
robust relative signatures under a frame shift since the PEMs have larger frequencies. 

The lowest order nonlinear terms come from the third variation of H expanded about the 
dynamical equilibrium, 
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In terms of the Fourier amplitudes this expression is, 

t 2 oo 

5 3 H = -Y V \—Mk k (p® aU) (h(J) - p U) aU) + c . c .) (28) 

j=l m,n=l 

+ ^7j(p {j) p {j) P ij) +P (j) p (j) p (j) +C.C.)]. 

Q A/f \' m 1 n ' —m — n 1 m 1 —n' n — m J J 

We then effect the same transformations on (28) that led to the diagonalized free energy (25). 
This spawns a myriad of nonlinear terms, only some of which survive an averaging process 
that leads to the Birkhoff normal form. 20 

For a three-wave resonance, one finds after near-identity transformations that the only 
higher order terms contributing to the normal form are of the type, 14 

0(3) Terms ~ J^ 2 jM/2 jM/2 gin ^ + m ^ + ^ ^ 

with |/| + \m\ + \n\ = 3. The Hamiltonian up to 3 rd order terms for the resonant NEM/PEM 
triplets can be written as, 



H = luiJ l — uj 2 J 2 — UJ3 J 3 + a^J J 1 J 2 J 3 sin (0 1 + 2 + 3 ) , (30) 

where a = a(k m , M, L) is a nonlinear coupling constant that is neither especially large or 
small in the parameter regime of the three-wave resonance considered here. 

Since the angles (9 1 ,9 2 ,9 3 ) appear in only one combination in H, further simplification 
of (30) is possible via the generating function, 

F 2 (1, 6) = X -l, {9, + 9 2 ) + I 2 6 2 + I 3 9 3 , (31) 

dF dF 

with ip =-*f- and J j =-q^- This canonical transformation yields, 

3 j 



a i 

H = uj, - u 3 I 3 + -y/l, (I, + 2/ 2 ) I 3 sin (2^ + V> 3 ) , (32) 

with 2u> 1 = lo 1 — u) 2 . If one chooses initial conditions satisfying I 2 = J 2 — J 1 = 0, then H is 
identical to the normal form of a two-wave interaction originally presented by Cherry. 21 ' 22 In 
terms of the (q, p) variables, Cherry's Hamiltonian is 

H = \^ (pI + Ql) ~ (f 3 + qfj + \ (2q lPl P 3 - q 3 [q* - $]) , (33) 
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where e = v ^ a . 

The dynamical system generated by the Cherry Hamiltonian is integrable. In the spe- 
cial case of a third-order resonance with uj 3 = 2u 1 , there exists a family of two-parameter 
solutions, 

q 1 = — —sin (0)^ + 77) , p 1 = -cos^t + r/) (34) 

g 3 = S in (20)^ + 27/) , p 3 = - cos (2^ t + 2r?) , 

where £ and 7/ are constants depending on the initial conditions. 

The solutions (34) show the possibility of finite-time density singularities when two neg- 
ative energy modes interact resonantly with a positive energy mode. A system exhibiting 
this behavior is said to undergo explosive growth, and it could be an important mechanism 
for structure formation in galactic and cosmological settings when relative motion between 
different fluid species is involved. If the resonance is detuned, separatrices bounding stable 
orbits emerge in phase space, but the dynamics are still prone to finite-amplitude instability. 

5 DISSIPATIVE INSTABILITY 

We close our investigation of the consequences of negative energy modes by examining 
the effects of dissipation on the linear stability of the system. With negative energy modes 
propagating through a dissipative medium, we may expect new instabilities since the damping 
can pump more negative energy into the wave. This somewhat counterintuitive effect of 
frictional forces in other contexts was first pointed out by Kelvin and Tait 23 (see also Zajac 24 ). 

Suppose that collisions are important at some stage in the development of a gravitation- 
ally bound structure. A simple model of this effect incorporates a dynamical friction term 
(— (m — 112) on the right hand side of the momentum equations (2), where v is a positive 
damping coefficient. 
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In the dimensionless symmetric case the dispersion relation becomes, 



1 - 2k 2 (M 2 + 1) lo' z + 2iv k\M A + 1) - 1 lo (35) 



.2 



.2/ ,,^2 



+k\M z -l) r(M 2 -l) + l =0. 

If we assume the damping is weak, !/ < 1, we may develop (35) in a regular perturbation 
series (o> = u; + + . . .) to find the lowest order corrections to the frequencies (10), 

u>t = -(l± 1 ^ , (36) 

2 V Vl-8fc 2 M 2 + 16M 2 W V ; 

where we assume k 2 » M± ^' 2 ~ - to avoid the singularity accompanying the vanishing de- 
nominator in (36) (see Casti 16 for the details). If we assume M > 1 so that negative energy 
modes are present, then a close examination of the corrections reveals that the dissipation 
promotes instability in the wavenumber band M± ^ 2 ~ - <C k < ^ for any M > 1 no matter 
how weak the damping. Since the instability as k 2 — > \ is realized only in the M — > oo limit 
of the undamped problem, we see that the dissipation indeed has the effect of destabilizing 
modes that were stable in the conservative case. A numerical investigation revealed that this 
result holds for any v > 0. 

The modal bands destabilized by the damping become more significant in the asymmetric 
case. As remarked in section 3.2, bubbles of unstable two-stream modes can pinch off from 
the Jeans- unstable bubble and result in well-separated instability bands. The inclusion of 
dissipation can destabilize the entire band of modes separating the bubbles, as well as some 
higher-fc modes beyond the undamped two-stream bubble. This is illustrated in figure 4. 

One should not assume that any form of dissipation will destabilize negative energy modes. 
For instance, if each fluid feels only a drag proportional to its own velocity, there are no new 
instabilities even with relative motion. In other words, the dissipation must in some sense 
project onto the eigenspace spanned by the negative energy modes in a way that decreases 
their energies. This depends not only on the nature of the dissipation, but also upon the 
initial equilibrium about which one perturbs. 
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The effect of damping can be understood by examining the time evolution of the free 
energy. If the dissipation acts to increase the energy of a positive energy mode or decrease 
the energy of a negative energy mode, then one can show that pure imaginary eigenvalues 
take on a positive real part. 25 To see that this is possible here, consider the temporal change 
in the free energy, which in the symmetric case can be written, 

5 2 H = -vM \ dx (5<p lx - 5<p2x) 2 -vM dx (Sfa - 6fc) (Sp lx + Sp 2x ) . (37) 
Jo Jo 

Since the first term of 5 2 H is negative definite, the conditions for which the free energy 
decays or grows in time is determined by the relative phasings of the velocity and density 
perturbations comprising the second term. One may deduce the effect of the dissipation on 
any particular mode of the conservative problem by inserting the undamped modes into (37), 
which yields a formula for S 2 H valid up to 0(y 2 ). For subsonic relative motion, M < 1, 
the expression (37) is negative definite and the damping lives up to its name and causes the 
PEMs to decay in time. When M > 1, 5 2 H can be either positive or negative for an NEM 
depending on the value of k, which explains why some NEMs are destabilized and others are 
damped in the usual sense. 

6 DISCUSSION 

The formation of structures through the action of gravity is much analyzed in cosmol- 
ogy, galactic structure and cosmogony. Most of this analysis is centered on the operation of 
gravitational instability, though streaming fluids can resonantly interact via the gravitational 
field to cause linear instability in spectral ranges inaccessible to the traditional Jeans insta- 
bility. As we have brought out here, the distinction between the two types of unstable modes, 
the Jeans and the two-stream, becomes sharper when one constituent is far denser than the 
other. Much of the previous work on the subject failed to take advantage of this crucial 
feature by focusing attention on situations where each component exists in equal abundance. 
Even when the two-stream instability does not occur, if the total energy of the gravitational 
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two-stream interaction is indefinite, the positive and negative energy modes that are stable in 
the linear theory can interact to produce explosive development of disturbances of arbitrarily 
small amplitude. This can be a significant aspect of the theory of structure formation. 

There are many clear instances where the dynamics of interpenetrating fluids may play a 
role in developing structures, but we close here by suggesting that even when the streaming 
is not apparent, two-stream dynamics may be relevant. An interesting example is provided 
by the coexistence of dark matter and luminous (baryonic) matter that is generally believed 
to occur throughout the cosmos. The locations of the two kinds of matter seem to be well 
correlated, which would not be the case if they were now streaming through each other. 
On the other hand, it might be reasonable to ask why there is this apparent correlation 
(or anticorrelation in the case of negative gravitational density) of the two kinds of material. 
Even if they had once been in relative motion this situation would not long persist, as we have 
seen. But the outcome, as far as large-scale structure is concerned, could be quite different 
if the kinematic history of the interaction of the two matters had been richer than has 
been supposed hitherto. Given the indefiniteness of the free energy if the initial streaming 
is large enough, waves of short length scale could have interacted in an explosive manner 
to quickly produce highly nonlinear density fluctuations. This is a feature of gravitational 
structure development that could be profitably studied, particularly in situations where the 
background Hubble expansion cannot be ignored. The dynamics of a two-fluid system with 
initially time- dependent relative motion is currently under investigation. 
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Symmetric Case , M = 2 
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Figure 1: Dispersion curves for supersonic motion, M = 2. The Jeans modes are indicated 
with a "J" and the two-stream modes with a "TS." The growth rates of the J and TS modes 
merge at the onset of oscillations shown in the lower panel. Dashed vertical lines demarcate 
the critical wavenumbers for growth of the Jeans and two-stream modes. The "+" and "— " 
signs in the frequency plot denote the energy signature discussed in section 4. 
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Asymmetric Case 




Figure 2: Asymmetric Case, (3 = .05, c = 1, M\ + M2 = 1.95. Illustration of the pinched 
two-stream bubbles. The frequencies corresponding to the two-stream modes are degenerate 
along the entire k— band of the bubble. At the edge of the bubble the frequencies once again 
separate. 
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Figure 3: Diagravic function in the symmetric case for linearly stable modes with k = 3. A 
stable mode exists wherever T crosses the w-axis. The positive energy modes are indicated 
with a "+" and the negative energy modes with a "— ". 
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Asymmetric Case With Damping 





Figure 4: Dispersion curves for damped supersonic motion, v = 0.1, M\ = 1.45, M<i = 0.5, in 
the asymmetric case with /3 = .05 and c = 1. The entire fc-band separating the Jeans bubble 
and two-stream bubbles of figure 2 are now unstable. The modal bands destabilized by the 
dissipation are demarcated by the arrows between the dashed vertical lines. 
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